Fano resonance in quadratic waveguide arrays 
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We study resonant light scattering in arrays of channel optical waveguides where tunable quadratic 
nonlinearity is introduced as nonlinear defects by periodic poling of single (or several) waveguides 
in the array. We describe novel features of wave scattering that can be observed in this structure 
and show that it is a good candidate for the first observation of Fano resonance in nonlinear optics. 



The study of nonlinear dynamics and spatial solitons 
in optical systems has recently attracted a great deal of 
attention y|. In particular, many specific properties of 
nonlinear lattice systems can be analyzed for arrays of 
weakly coupled optical waveguides where both nonlin- 
earity and diffraction may differ dramatically compared 
to those in the corresponding continuous systems 0, El • 

During last years a growing interest is observed in the 
study of nonlinear optics associated with the so-called 
quadratic nonlinearities which may produce the effects 
resembling those known to occur in cubic nonlinear mate- 
rials. Typical examples are all-optical switching phenom- 
ena in interferometric or coupler configurations as well as 
the formation of spatial and temporal solitons in planar 
waveguides 4]. Recently, it was demonstrated experi- 
mentally PI that arrays of coupled channel waveguides 
fabricated in a periodically poled Lithium Niobate slab 
represent a convenient system to verify experimentally 
many theoretical predictions, including the first observa- 
tion of two- frequency discrete solitons mutually locked by 
quadratic nonlinearity. These experimental observations 
open many perspectives for employing larger nonlineari- 
ties in lattice systems made of quadratic materials. 

In this Letter we suggest to employ the arrays of 
weakly coupled nonlinear quadratic waveguides for the 
study of novel effects in resonant light scattering. In par- 
ticular, we show that when periodic poling is applied to 
just a few waveguides in the array, it creates a nonlinear 
defect [H that possesses specific resonant scattering 
properties and may be employed for the first experimen- 
tal observation of Fano resonance in nonlinear optics. 

Following the waveguide design recently implemented 
in Ref. Q for the observation of discrete quadratic soli- 
tons, we consider a discrete model describing an array 
of weakly coupled linear waveguides where one or sev- 
eral neighboring waveguides have periodic poling and 
therefore possess a quadratic nonlinear response (see 
Fig-O- When the matching conditions are satisfied, the 
fundamental-frequency (FF) mode with the frequency u 
can generate paramctrically the second-harmonic (SH) 
wave at the frequency 2u> , so that such a structure with 
several poled waveguides may behave as a nonlinear de- 
fect with localized quadratic nonlinearity. The continu- 



FIG. 1: Schematic of an array of channel waveguides with 
quadratic nonlinearity where nonlinear defects are created by 
periodic poling. Arrows show the scattering process. 



ous version of this model has been studied earlier [8j . 

In the tight-binding approximation usually employed 
in the theory of discrete lattices 0, the effective equa- 
tions for the complex envelopes of the FF wave (u n ) and 
its SH component (v n ) coupled at the defect waveguides 
with n — 0, . . . , N can be written in the dimensionlcss 
form 



,du n 



N 



i-p 1 + c u (u n+ i + U n -i) + 2 ^2 u * m Vm5 nm = 0, 



dz 
. dv. 



N 



(1) 



i-^j 1 + Cv{v n +i + v n -i) - Av n + ^2 u 2 m 5 nm = 0, 



where c u and c v are the coupling coefficients, S nm is the 
Kronecker symbol and A is the phase mismatch param- 
eter assumed to be identical for all waveguides. 

First, we analyze the scattering of a plane FF wave 
at the frequency lu by the a single (n = 0) quadratic 
defect waveguide ('impurity site'). After the interaction 
with the quadratic waveguide, the FF wave generates a 
SH wave which could either propagate or get trapped 
being guided by the defect waveguide. To calculate the 
transmission coefficient t(k) of the FF wave, we present 
the fields in the form, 



u n (z) = exp(i/3iz) 



J /exp(ifen) + i?exp(- 
1 Texp(ifcn) 



ikn) n < 
n > N, 
(2) 
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v n (z) = exp(i/3 2 z) 



J Rexp(-iqn) 
\ Texp(iqn) 



n<0 
n> N, 



(3) 



where (3% = 2c u cos k and /3 2 = 2c„ cos g — A are propaga- 
tion constants of the FF and SH respectively, and k and 
q are corresponding transverse wavenumbers. By using 
the phase-matching condition Q {2fi\ = /?2) we obtain 
the relation c v cosq — (A/2) = 2c u cosfc, that defines the 
dependence q = q(k). For k m - m < k < fc max , the function 
q(k) takes real and positive values. Outside this interval, 
the values of q(k) are purely imaginary, and they corre- 
spond to localized (non-radiating) states trapped by the 
defect waveguide at n = 0. A simple calculation yields 
the result: &: m ax,min = cos _1 [(c„ ± A/2)/2c„], when these 
values are real and positive or zero, otherwise. 

Evaluating the mode coupling at the impurity (?i = 
0) and neighboring (n = —1) sites allows to obtain 
the relations, T = I + B and R = T, and derive a 
nonlinear equation for the transmission coefficient t = 
\T\ 2 /\I\ 2 in the form, t[l + b(k)t} 2 = 1, where b(k) = 
\I\ 2 [2c u c v sinfcsing(fc)] , that has only one real solu- 
tion. The resonant scattering, when a localized SH field is 
generated, can be analyzed similarly by replacing q — * iq. 

The study of wave scattering in this system predicts 
the resonant suppression of transmission at some points, 
i.e. i(A minjmax ) = (see Fig. EJ. We demonstrate below 
that these resonant reflections correspond to a novel type 
of the well-known Fano resonance Indeed, accord- 
ing to the Fano theory [9|, destructive interference and 
resonant suppression of transmission are observed when 
there exists a localized state coupled to the propagating 
channel with the energy inside the linear spectrum. Note, 
that q(fc m i n ) = and g(/c m ax) = 7r, i.e. these values de- 
fine the band edges of the propagation spectrum of the 
propagating SH field, and the resonances take place when 
the SH field is generated. This situation seems to be in 
contradiction with the classical definition of the Fano res- 
onance. However, below we demonstrate in more details 
that this kind of the resonant scattering can be indeed 
defined as being associated with the Fano resonance. 

First, we consider the simplest case when the coupling 
between the SH modes of different waveguides vanishes, 
i.e. c v = 0, which is indeed the case of the recent ex- 
periments [f|. Then, in the stationary case, the coupled 
equations can be written in the form 



(3iu n = c u {u n+ i + u n -i) + 2uQv8 n0 , 



2f3iv = -Av + ul, 



(4) 



This model Q describes the main propagation channel 
for the field u n and an additional discrete mode v = v a 
coupled to it parametrically, and it is similar to the so- 
called Fano- Anderson model 01 1 except that the cou- 
pling here is nonlinear, that makes the scattering prob- 
lem nonlinear. 

To simplify the analysis, we eliminate the discrete 




FIG. 2: Examples of the FF transmission coefficient for the 
scattering by the quadratic impurity waveguide, for c u = 1, 
and (a) c v = 0, A = 0, (b) c v = 0, A = 2, (c) c v = 0.5, A = 0, 
and (d) c v = 0.5, A = 2, for two values of the intensity: 1 = 1 
(solid line) and 1 = 2 (dashed lines). 



mode described by the second equation and obtain, 
a i i \ , 2K| 2 u ; . 

PlU n = C u (U n+1 + U n -l) + — — -OnO, 

Api + A 



(5) 



which is an effective equation for the propagation channel 
that contains a scattering potential. The strength of this 
nonlinear resonant scattering potential depends on the 
incoming intensity \I\ 2 and propagation constant /3i(fc). 
If A is chosen such that kp is between and tt, and 



2Pi(kt 



A, 



(6) 



our potential becomes infinitely large for a particular fre- 
quency Pi(kp), which will lead to the perfect reflection. 



Note here, that kp 



v ffliii,max 



for our case c„ 



0. In- 



deed, after some algebra we can write down the equation 
for the transmission coefficient in the following form 



t 3 + 7 W(t-l) = 0, 



(7) 



where j(k) — (4c„ cos k + A) 2 c 2 sin 2 k/\I\ 4 . From this 
equation one can see that, when j(k) — then transmis- 
sion coefficient t vanishes. This happens at the wavenum- 
bers k = and k = tt, which correspond to the band 
edges of the propagation spectrum, and also at k = kp. 
In the latter case, this is exactly the Fano resonance. 

When the coupling between the SH modes in the 
waveguide array does not vanish (i.e. c v ^ 0), it leads 
to the appearance of the spectrum of propagating SH 
modes, /?2 [<?(&)]. At the band edges of this spectrum, 
/?2(0) and /^(tt), which correspond to k = k nlln and 
k = k max , the propagating SH field is described as a 
standing constant-amplitude mode of the forms v n = vq 
and v n = (— l) n VQ, respectively, where vq is constant. 
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FIG. 3: Top: Resonant generation of the trapped second- 
harmonic field by a Gaussian beam of the FF (a) and SH (b) 
fields for the parameters c u = 1, c„ = 0, A = 0, and k — 
tt/2. Bottom: Comparison of the transmission coefficients 
of the plane waves (solid) and Gaussian beam (crosses) for 
c u = 1, c„ = and (c) A = and (d) k — tt/2. Transmission 
coefficient of the Gaussian beam does not vanish due to its 
finite spectral width. 

Therefore, for these two cases we can again obtain the 
single-site equation for the second scattering channel. By- 
applying similar approach to these particular cases, we 
obtain two conditions for the Fano resonance, 

2/3i(fc Flia ) = -A±2c, (8) 

which occur exactly when the propagation constant 2/3 of 
the generated SH field coincides with either the propaga- 
tion constant /^(O) or fain), i.e. at the band edges of the 
linear spectrum of the propagating SH modes. In other 
words, in such situations we excite 'constant modes' by 
a local perturbation. Since the group velocity of these 
modes vanishes at the band edges, any local excitation 



could not propagate at the given frequency. It makes 
these modes effectively local and leads, finally, to the phe- 
nomenon of Fano resonance. 

We note here that, in our physical system of quadratic 
waveguides, the coupling between the first and second 
propagation channels is nonlinear and, therefore, it de- 
pends on the intensity of the incoming wave. According 
to the formulas © and (JSJ the position of Fano reso- 
nances does not depend on the value of this coupling due 
to its local nature As a consequence, this novel type 
of Fano resonance should exist for any intensity of in- 
coming waves similar to the conventional Fano resonance 
in the linear theory. But the width of the resonance de- 
pends on this coupling|l2j and, therefore, on incoming 
intensity of light (see Fig. 

In order to check the validity of our plane wave analysis 
and the manifestation of the effect in a realistic experi- 
ment, we have performed the numerical simulation of the 
Gaussian beam scattering. The results are summarized 
in Fig. |3 and they are in a good agreement with the 
theory of plane wave scattering. 

In the case of N defects and vanishing coupling be- 
tween them (c v = 0), Fano resonance does not change 
its position © [01 . After scattering by the first defect 
close to Fano resonance, other waveguides become almost 
transparent due to a small incoming intensity. Therefore, 
the width of the resonance will remain almost the same 
as in the case of a single defect. 

In conclusion, we have analyzed a novel waveguide 
structure where an optical analog of the Fano resonance 
can be observed as peculiarities of the resonant scattering 
and second-harmonic generation in the quadratic waveg- 
uide arrays. We believe this kind of waveguide struc- 
tures, already fabricated for the observation of discrete 
quadratic solitons ||, is a good candidate for the first 
experimental observation of Fano resonances in nonlin- 
ear optics. 
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